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Abstract 

We study symmetries of the Bethe equations for the Gaudin model appeared naturaUy in the 
framework of the geometric Langlands correspondence under the name of Hecke operators and 
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Introduction 



The Gaudin model is one of the first nontrivial examples of quantum integrable systems which shows 
deep connection of problems of mathematical physics and such modern fundamental mathematical 
questions as the Langlands correspondence. For CP^ with marked points the geometric Langlands 
correspondence is an interpretation of the separation of variables in quantum Gaudin model. The 
basic purpose of this work is to specify the Langlands correspondence concerning the problem 
of solution of a quantum model. We use the form of the Bethe equations which expresses a 
condition for the spectrum of the Gaudin model in terms of the monodromy properties of certain 
differential operator with regular singularities (called the G-oper). Further we construct the discrete 
monodromy preserving transformations known as Schlesinger or Hecke transformations. It appears, 
that they determine also an action on rational solutions of the G-oper, and therefore on eigenvectors 
of the Gaudin model. It is necessary to mention that such transformations act by changing the 
highest weights in representations and also the inhomogeneity parameters of the model. 

The exposition is organized as follows: in section [T] the Gaudin model is defined, the traditional 
version of the Bethe ansatz and its formulation in terms of monodromy properties of specific 
differential equations is described. In section [2] the "matrix" form of the monodromy condition 
equivalent to the Bethe equations is constructed. In section [3] the Schlesinger transformations are 
investigated. The last section is devoted to some relations of this subject with other problems. 
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1 Bethe ansatz for the Gaudin model 
1.1 Gaudin model 

The Lax operator of this model is a rational function with simple poles: 



where E^i form the standard basis in Main and e^/ form a basis in the i-th copy of the Lie algebra 
g[„. It is the same as to say that the quantum algebra is 




N 



with the residues G Main <8) fltn C Mat„ ® ?7(0[„)'^^ defined by the formula: 




(1) 



kl 



and the commutation relations for the Lax operator are 



[L(z) 1,1® Liu)] = [Ri2{z - u),L(z) (g>l + l(g) L{u)] G Mat^^ (g> A 
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Let us now define the "quantum" determinant of a matrix with non commutative entries as the 
completely symmetrized determinant 

"det"(i?) = (-l)'^'^-Sr(l),<T(l) • • • B^(^n),a{n) 

T,(TgE„ 

Then the quantum spectral curve is defined as follows 

n 

"(iet"(L(z) - d,) =^Qh{z)dr' (2) 

k=0 

Theorem [T] The coefficients QIk{z) commute 

[Q/fe(z),Q/^(n)] =0 
and quantize the classical Gaudin hamiltonians. 

Remark 1 By the same formula in |^ it was constructed a commutative subalgebra in C/(0l„)[t]/i^, 
J7(g[„)[t] and the center of the universal enveloping affine algebra on the critical level . 

1.2 Bethe ansatz 
1.2.1 sl2 case 

Let us consider the Gaudin model for the 5I2 case. The Lax operator can be expressed by the 
formula 

N 



_ I A{z) B{z) \_sr^^ 
\ C{z) D{z) j ^^z-z,. 



i=l 



where 

^ _ , hi/2 Ci 



f^ -hi/2 

Let us also fix a finite dimensional representation V\ = ®Vx^ which is the tensor product of 
highest weight representations for each factor f7(5(2) of the quantum algebra. The vacuum vector 
is \vac >= ®j|0 >i where the vectors |0 are defined by ej|0 >i= and hi\Q >= AjjO > . The 
quantum characteristic polynomial in this case is 

"def'{L{z) - d,) = -^TrL\z) + 9^ 



TrL^(z) - V ^"^^ + 2 V ^"(^-^^-^ 



Kj 

The coefficients at the double poles are central 

rr$2 = + CiJi + fie, = 4 
and take values l/2Ai(Aj + 2) in representation V\. Hence 



2 ^ (z - Zi)^ ^ z - Zi 
where the residues Hi are the Gaudin hamiltonians 

H - ^ ^»^i/2 + e^/j + ejfi 

. , . Zi Zj 
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1.2.2 Bethe vector 

The Bethe ansatz method supposes that a common eigenvector Q in V\ 

HiQ = xnmn (3) 

has the form 

M 

Q = Y\_C{fij)\vac > 

i=i 

for some values of parameters . This condition imphes the following system of algebraic equation 
on parameters iJ,j (the Bethe system) 

-^E^^ + E^— = ^■ = i'--->^ (4) 

The eigenvalues take the form 

A 



\ 3 ■' J7^« 



•3 



Zi Zj 



1.3 "Null-monodromic" form 

Let 

^{z) = \{{z-z,r^^'^J\{z-^,,) (5) 

be a solution of the equation 

/ 1 „(2) rr \ 

"def^iL{z) - d.Mz) = l^d', - 2 E (731^ - E ^ j ^(^) = (6) 
supposing that Hi are scalars. This is equivalent to the following system of equations 



1 , / V- 



Z Zj \ Zj Z T Zj LL j 

1 _ 1 \i ^ 1 







which is exactly the system of Bethe equations and the conditions on eigenvalues. Remarkably the 
second solution of the equation [6] has the same behavior at z = Zi, that is at most half integer 
exponents. This imply that the monodromy of the fundamental solution lies in Z/2Z. There is also 
an inverse statement 
Theorem [4] The equation 

with z/j G C has monodromy ±1 iff i^i are equal to xn(-f^i) for some common eigenvector Q of the 
Gaudin hamiltonians in the representation V\. 
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2 Matrix form 



2.1 Scalar versus matrix G-opers 

Let us now consider a connection in the trivial rank 2 bundle on the punctured disc of the form: 



k 



aii{z) ai2{z) \ _sr^ A 
a2iiz) 022(2 

where the residues satisfy the following conditions 



anz «i2. ]=y^ (8) 
V a2i{z) 022(2) J ^ z-Zi 

1 = 1 



Tr{A,)=0; Det{A,) = -dh E^^=(o -k) ' 

Let us now consider the Fuchsian system 

{d,-A{z))^{z)=0 (10) 

which is equivalent to the following 

ip[ = auTpi + ai2Tp2 

-02 = a2lV'l + ^22^2 

This system implies 

V'l = a'i2/ai2'>P[ + uipi 

where 

u = a[i + aji - ana[2/ai2 + 012021 
Then if one passes to $ = V'l/x where x = \/ai2 one obtains the equation 

+ ;7^> = 

where 

= x'Vx - a'i2/ai2x7x - (H) 

More precisely this implies 

U=U^)--(^)\an^- a'n - o?, - 012021 (12) 
2 V«i2/ 4 V«i2/ ai2 

Suppose that 012(2) has distinct zeros 



nk-2f 
j=l[Z-Wj 

ai2{z) = c—jr 



Ui=iiz-Zi) 

Let us note that the number of zeros is in agreement with the normalization condition [H Hence 

k-2 , k 



^ = y^ — y^— (13) 



012 ^ Z-Wj ^ 2 - Z, 
j = l 1=1 
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The expression for the potential takes the following form 

^(Z-Wj)^ ^ {z-ZiY ^Z-Wj ^Z-Zi 



where 



Hi, 



1 i \ V- 1 V- Tr{AiAj) + all + <\ + 1/2 
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Let us note that the coefficients at (z — Zj) ^ can be rewritten in the form 

I / 4. + detAi = {1/2 -di){l/2 + di) (15) 
Further we will identify this with the values of Casimir elements. In this way Aj = 2di — 1. 

2.2 Dual G-oper 

As it was shown in previous calculation the matrix connection is related to the Sturm-Liouville 
operator with additional regular singularities at points Wj. The similar consideration for the second 
component ^'2 can produce another scalar differential operator of the second order with poles at Zi 
and additional points Wj defined by the formula 



021(2) 



]ti=i{z-Zi) 

We will call the corresponding Sturm-Liouville operator 

dl-U 

the dual G-oper where 
2.3 Pull back 

In this section we construct the inverse map, namely from the scalar G-oper, which has trivial 
monodromy in the sense of the Bethe equations, we construct a rank 2 connection of the form 1101 
such that its horisontal sections have monodromy in Z/2Z. 

Let us consider the solution for the matrix linear equation [TOl 

{d,-A{z))^ = 

in the following ansatz 

k 

i^l = ll{z-z.,)-'^Mz) (17) 

i=l 
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where 

M 
M 

<^2M = ^^^ (18) 

Let us rewrite the system [10] with respect to the new parameterization [TSl 

dztpi/tpi = au + ai2</'2M (19) 

dzi^l/lplih/h) + dz{4>2/4>l) = 0.21 + a224>2/4>l (20) 

Let us represent these equations in a more explicit form: 

_ ^» _^ 1 ^ "-11 _^ °i2 y^ ^21) 

^-^ z — Zi ^-^ z — 7,- ^ z — Zi ^ z — Zi ^ z — 7,- 

i j i i- j 

y^ '^^ _^ y^ ^ j y^ "i _ y^ ^ y^ "21 _ y^ "n y^ "j ^22) 

2 — ^-^ z — I ^ z — 7,- ^ (z — 7,-)^ z — Zi ^ z — Zi ^ z — 7,- 

Further we find the conditions, that the residues of both sides of I19|20l at z = Zi coincide: 

11 + «l2 E (23) 



a' 



-Et^^^ = «--«nE7^ (24) 

These equations together with the trace condition + = imply that Si are the eigenvalues 
of Ai, in particular the choice Si = di is appropriate. Let us consider the behavior at poles z = 7j- 
Remarkably the poles of the second order at these points in equation [23] reduce. Let us calculate 
teh residues of both sides of the equations [23] and [2l] 



y-^^ + V^_ = (26) 

fe^^'-^V fe^^'"^^ i ^^--^^ 

Let us recall the normalizing condition on the residue at 00 

k 

E«i2 = (27) 
1=1 
k 



E«2i = (28) 



■t=i 



Let us note that the choice of poles for the Sturm-Liouville operator fixes the zeros of the rational 
function ai2{z), which in turn is defined up to a constant: 



012(2;) = c 



Y[)=1{Z - Wj) 

Ui=iiz - Zi) 
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Then the condition [27] fulfills automatically. The coefficients a\2 can be expressed by the formula 

nj=i(^i - 

The coefficients a\i can also be expressed as follows in virtue of 



a\, = c'-^^ ^ (29) 



rfc-2 

1. 

— ~Si — C 



Uj=iizi-w,) ai 



Let us substitute the expressions for a\2 and in equations [25l [26l and express then from the 
first one and substitute to the other: 

Wjin-j - Zj) I _ 2Sk ^ 1 ^ Uijzk - Wl)l[sM7ni - Zs) 

Ui{7j-Wi)y ^Ij-Zk ^.Ij-lk l-;^Y\l^k{zk- Zl)Y\s{lra-Ws){lj- Zk) 

Let us pass to the equivalent form dividing by yi'^^^^-w\ 

^ 2sk 1 ^ Uiizk-wi)UsM^'^ 



k ij-zk ^.ij-ik f^n/^fc(^fe-^i)n.(7m-«^s)(7i-^fc) 

n»(7fc - Zj) Wmilj - Wm) _ Q 

^ Yiiilk - Wi){-fj - -fk) Umi^j - Zm) 

Consider the left hand side of the equality as a rational function F{^j) and calculate its simple 
fraction decomposition at poles z^, w^, "fk and oo. It turns out that this decomposition takes the 
form: 

Y7i-^fc ^7J-Wk t^j^j-^i 

Hence, the explored equality is equivalent to one of the Bethe equation. Let us demonstrate, for 
example, the calculation of the residue at 7j = Wi 

Res Fh-) - V nlizk-wOUsM^^-Zs) 
Wsiwi-Zs) liiizk-wi) 

Mm: 



Wsi^ik^i -'^■■>)k n/^fc(^fc - Zl)iwi - Zk) 

Let us note that the expression on the right is of the form 

{Res,=^Mz)r^Y.^^'^=^>'^^ 

k 

where 

11/ - wi) 



Hz) 



Wliz - Zi){z - Wif 

and hence is equal to —1. 

Then the sufficient condition expresses in the following 
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Theorem 1 If the set of numbers 7j where i = 1, . . . , M satisfy the system of Bethe equations with 
parameters: the set of poles is zi, . . . , z/^ and wi, . . . , with the highest weghts 2si — 1, . . . , 2sfc — 1 
and correspondingly, then the vector 

i=l 

where 

M 



h/<t>i = Y.^- (34) 
and the coefficients aj given by the expressions 



(35) 



solves the matrix linear problem [721 where the connection coefficients are defined by 

nAZi — Wn) 

= n,..(^. - '''' 

the coefficients a\i and a^i are defined from The normalizing conditions O hold. 

Proof. Properly, we need to that the normahzation condition [28] does not depend on the choice 
of the parameter c, in particular can be taken to be equal 1. Indeed, in virtue of 1231 [2^ we obtain: 

„.,_L5:^+i!4^fy:^V'i (37) 

We need to prove that 

The first summand of [38] can be transformed to the following expression using the Bethe equations: 

^ '1^ ^l3-'^ ^1,-^^ fe^^--^V 

Let us simplify the second summand of [38] changing the summation order 



(40) 
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Considering the second summand of[30]let us emphasize that 



where 



j^i {Zi -Zj){Zi--im) n j (7m - Zj ) 

Hence the expression 1411 takes the form: 



(42) 




WAlm - Zj) \^ 7m - W. 



(43) 



which reduces the corresponding part of 1391 Let us consider the first summand of 1401 it can also 
be simplified: 

Wji^i-Wj) _ \[j{ll-Wj) Wjilm-Wj) 

^ Wj^iiZi - Zj){Zi - -im){zi - 7/) Y\j{ll - Zj){jm - ll) Wji^m - Zj){jm - 7/) 

Inserting this expression in |30] we finish the proof □ 

3 Schlesinger transformations 

There is a discrete group of transformations which preserve the form of the connection II 1 and more- 
over do not change the class of monodromy representation. However they change the characteristic 
exponents at singular points by half-integers in a specific manner. Such transformations are called 
Schlesinger, Hecke or Baklund transformations in different contexts. This type of transformation 
has simple geometric interpretation. 

3.1 Action on bundles 

Let us first of all consider holomorphic bundles on CP^. Due to the Birkhof-Grotendieck theorem 
any holomorphic bundle on CP^ of rank n is isomorphic to 0{ki) © ... © 0{kn) for some tuple 
{ki, . . . , kn) which is called the type of the bundle and is defined up to the symmetric group action. 

Let us now consider a covering of CP^ which consist of C/qo - a disc around oo not containing 
z = Zi, i = 1, . . . , N and Uq = CP^\{oo} . We will work with rk = 2 holomorphic bundles and 
parameterize them by the gluing function G{z) - holomorphic invertible function in Uq D C/qo with 
values in GL{2). We say that the pair Soo{z) G O^^^Uoo) and Sq{z) £ 0^'^\Uo) defines a global 
section if So{z) = G{z)Soo{z). 

All transformations will be described by some change of the gluing function. Let us consider 
such a transformation which multiplies the gluing function on the left by 

Gs{z) = Gs(^'~'' (45) 
for some constant matrix Gs and some point Zs G Uq. 
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Remark 2 The action on gluing functions can be brought to the action on isomorphism classes of 
bundles making the choice of constant matrix Gg dependent on the trivialization, more precisely if 
one change the trivialization in Uq by T(z) one should change the matrix Gg by T{zs)Gs. This is 
clear from the remark\^ on the invariant definition. 

The case of our interest will be the composition of two such transformations at different points. 



Lemma 1 The composition of two transformations given by Gi{z)Gj^{z) applied to the trivial 
bundle give also the trivial one for generic choice of constant matrices Gi, Gj. 

Proof Our goal is to perform the decomposition G{z) = Gi{z)Gj^{z) = Gij{z)Goo{z) where 
Gij{z), Gooiz) are invertible respectively in C/q, Uqo. The general argument here is the observation 
about the dimension of cohomology spaces in families at generic point. Indeed, for particular choice 
of G~^Gj = 1 one obtains the trivial bundle which is semistable and hence minimizes the dimension 
of H^[End{V)) for V of degree which is the case. Due to this the trivial bundle is generic in the 
family of bundles given by different G. 

However we provide here an explicit demonstration in spirit of the decomposition lemma of [5]. 
Let us decompose the product 



Giz) = G^^^ ^ j GT'G, ^ ^ Q ^ ^ J GJ^ (46) 
where we use the following notations 

= ( ^. " ) '''' 

into the product 

G{z) = G^j{z)G^{z) 

where Gij{z), Goo{z) are holomorphic invertible functions in Uq, Uoo respectively. Traditional 
calculations provide a decomposition of the form 



Gooiz) 



z{l-Xjyi){l-Xjyj)-Xj(y^-2yj-Xjyiyj) Xj 

(l-2xjyi+Xjyj){l-Xjyi){l-Xjyj){z-l) {l-Xjyi){l-Xjyj){z-l) 

yi-2yj+Xjyiyj 1 

{l-Xjyi){l-2xj yi+xjyj) l-Xjyt 



3.2 Transformations on connections 

Let us define an action of the group of Schlesinger transformations Tij on the space of connections 
with singularities at z = Zi. As mentioned the action on bundles changes the gluing function. Let 
us start with the trivial bundle given by the gluing function 1 . The Schlesinger transformation gives 
another bundle structure, a global section is defined by a pair 5o,5oo such that 5*0 = GSoo with 
G = GijGoQ- One can define an action on the space of connections as follows: let dz — A he an 
initial connection in the trivial bundle defined by the same expression over two considered open 
sets, the transformed connection is then defined as follows: 

dz — A over Uoo 
G{dz-A)G-^ over Uq 

After changing basis in Uoo of the form Soo = GooSoo we obtain the expression for the connection 

dz-A ^ Goo{dz-A)G^' (48) 
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If one changes the triviahzation in Uq as follows 5*0 = G-^Sq one obtains 

G{d, - A)G-' ^ GT^'G{d, - A)G~'Gi, = Goo(9, - A)G^' (49) 

Hence the transformed connection can be represented as a global rational connection of the same 
type as the initial one, the same behavior at oo is guarantied by the fact that Goo is holomorphic 
invertible at Uoo- 

Remark 3 This definition is coherent with the invariant construction, the Hecke transformation 
on the sheaf language is the operation of taking a subsheaf, defined by a linear condition on values 
of sections at fixed point. To define an action on the space of singular connections at a singular 
point one has to choose as a condition the dual covector to one of the eigen direction for the residue 
of the connection. 

Using the results of the previous section let us calculate the action of the Schlesinger trans- 
formation on the connection. To observe the normalization condition of A{z) at oo one should 
consider the transformation of the form 



G{z) = G^i(oo)Goo(^) 



1 



1 



XI (yo - 2yi+x i ypyi ) xi(l~23;iyo+a:ii/i ) \ \ 

{l-xiyo){l-xiyi) (l-xiyo){l-xiyi) \ \ (r.r^\ 

yo-2yi+xiyoyi l-2xii;o+zi^/i I I V"-"-"/ 

{l-xiyo){l-xiyi) (l-xiyo){l-xiyi) / / 



The whole set of Schlesinger transformations for the 3-point case relevant to the analysis of the 
Painleve VI equation was calculated in [6| 

3.3 Action on Bethe vectors 

Let us introduce the notation Z = {zi}, W = {wj}, A = {Xi}. Let us also introduce symbols for 
the following objects: 

• OS{Z, W,A) - a solution for the system of Bethe equations, that is the set fii which solves 
the system H] for the fixed points (the poles of the Lax operator) zi, . . . ,Zk and wi, . . . , Wk~2 
and the corresponding highest weights Ai, . . . , Afc and 1, . . . , 1; 

• OL{Z, 1^, A) - a solution of the linear problemOof the formOwith the given set of parameters; 

• ML{Z,W, A) - a solution of the matrix linear problem 1101 of the form [T7] with the given set 
of parameters, it has two components ML{Z, W, A)i^2] 

Remark 4 The choice of the highest weights 1 at "moving poles" Wi is not necessary but in a sense 
generic. One can take the potential 

U = V -'/"'"^ +f + f + "-y- + g + f -gf^ (51) 

'f—f (Z — Wj)''' ^ [Z — ZiY ^ Z — Wn ^ Z — Zi 

J=l ^ J' 1=1 ^ '■' J=l -I 1=1 ' 

with generic values of the highest weights . This realizes if one demand ai2iz) to have zeros Wj 
with multiplicities r]j subject to the condition X^"!]^ r]j = k — 2. 



The results of our previous considerations can be summarized in the following: for a solution of 
the Bethe system OS{Z,W, A) one can consider a solution ML{Z,W, A) in virtue of the theorem 
[TJ Then one can define a family of transformations: 
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• Dual solution Symbolically this transformation expresses vie the diagram 




• Schlesinger transformation Using the notation below one has 




Depending on the choice of the subspace of the upper and lower Schlesinger transformations 
one obtains the following action of Tij on the set of highest weights 



Discussion 

Painleve equations 

The idea of discreet transformations preserving the condition of "minimality" of the monodromy 
exploited here lies in the middle of the theory of isomonodromic transformations. The same group of 
transformations act on the space of solutions of the Painleve VI equation realized as the monodromy 
preserving condition for the Fuchsian system (the references can be found for example in [7]): 



One can expect that the problem of continuous deformation of the spectrum for the quantum 
Gaudin model with respect to the parameters of the model - the marked points, is equivalent to the 
Schlesinger system of equations, hence for some particular case of the model is just the Painleve 
VI equation. This is not a priori evident because of the essentially resonant case of the connection 
appearing in the isomonodromic formulation of the Bethe equation for the Gaudin model. 

Langlands correspondence 

The quantum Gaudin model plays an important role in the geometric Langlands correspondence 
[8l[3], namely the correspondence is provided by the relation: 



In the arithmetic case moreover one has an important property of coincidence of the corre- 
sponding L-functions. In the geometric case over C there is no natural notion of the Probenius 
automorphism and hence of the Galois-side L-function. Presumably, the symmetry elaborated here 
is related namely to the Galois side of the correspondence. 

Recently in [9] it was observed another appearence of the arithmetic L-functions in quantum 
integrable systems, namely in the theory of quantum affine Toda chain. 




(9, - A{z))^{z) = 



where 




Gaudin-Hitchin" D-module 4^ flat connection (G-oper) 
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